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Abstract 

When a vector meson such as J/'E'(35) or Y(4S') decays to a particle- 
antiparticle pair of neutral mesons, the time distribution of inclusive decay 
to a given final state is naively expected to be the incoherent sum of those of 
the two mesons with opposite flavors. In this paper, we show that this is in 
general not the case for arbitrary coherent two-body states of the mesons, and 
obtain conditions under which such a naive incoherent sum gives the correct 
distributions. The analysis is based on the Weisskopf-Wigner formalism, and 
applicable to the cases where there are more than two orthogonal states that 

can mix to form a set of eigenstates of mass and decay rate. 
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I. INTRODUCTION 


In the studies of CP violation in neutral meson systems such as 0 —> K°K° [!1|, 
J/'h(3S') —> D°D°, or T(4S') —> B°B° one often deals with an inclusive decay dis¬ 
tribution where one of the mesons decay to a given final state / at a given time t and the 
other meson can decay to any final state at any time. It has been shown in a recent study 
that in the decay of T(4S'), the inclusive decay time distribution of T(45) —► / is the in¬ 
coherent sum of the decay time distribution of a pure B° at t = 0 decaying to / at time t 
(denoted T B o^f(t)) and that of B° (denoted V B o^f(t)) 0. Such relation is critical in anal¬ 
yses of inclusive lepton asymmetries [[3| , or in relating observed inclusive yield asymmetries 
of Y(4S') —> / and T(4S') —> f to the asymmetry of decay amplitudes Amp(B° —> /) and 
Amp(B° —> /) H]. When a T(4S') decays to a pair of neutral B mesons, it is generated in 
a coherent L = 1 state, which is antisymmetric under the exchange of the two mesons: 

4= (|B°M}|B 0 (-h> - |B°(b>|B°(-fc)}) . (1-1) 

where the mesons are labeled by their momentum ±k which we will drop hereafter and 
implicitly assume that left (right) side of the meson pair is in +k(—k) direction. A natural 
question is then whether such an incoherent sum gives the correct inclusive distribution for 
a general two-body state given by 

$ = aB°B° + bB°B° + cB°B° + dB°B° , (1.2) 


where a, b, c, d are arbitrary complex coefficients with 


\a\ 2 + \b\ 2 + | c| 2 + \d\ r “ 


1 . 


When the pair is generated with a definite orbital angular momentum, further symmetry 
relations must be satisfied; in this study, however, we will keep the general form as above. 
As we will show below, the necessary and sufficient condition for the naive incoherent sum 


(2|a| 2 + \b\ 2 + |c| 2 ) T B o^ f (t) 

+ (2|d| 2 + |6| 2 + |c| 2 ) r B o^ f (t) 


(1.3) 
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to give the correct distribution for any final state / (and independent of the details of the 
mixing) is 

G = a*(b + c) + d(b + c)* = 0. (1.4) 

II. GENERAL COHERENT TWO-BODY STATES 

In the following, we study a system of n orthogonal states Bi (i = 1 mixing 

to form n eigenstates of mass and decay rate (physical states) B a (a = a,b,c...). The 
eigenstates B a are not necessarily orthogonal when CP is violated. We use the Weisskopf- 
Wigner formalism ||, but no assumptions are made on CP or CPT symmetries unless 
otherwise stated. The essential approximation used in the formalism is that the oscillations 
caused by mass differences and the decay rates are sufficiently slower than the time scale of 
decay transitions, which is a very good assumption for the cases under study 0 . 

The eigenstates B a then evolve as 

B a - e a (t)B a , e a (t) = e -^/ 2 +-W ? (2.1) 

where 7 a and m a are the decay rate and mass of the physical state B a . The time t is the 
proper time of the particle under consideration. The Weisskopf-Wigner formalism can be 
relativistically extended to moving particles; it can be shown, however, that it is equivalent 
to the evolution in the rest frame formulated as above 0 . We will hereafter consistently use 
the indices i,j for the orthogonal states B±, B 2 ..., and Greek indexes a, (3 for the physical 
states B a , B b ... : 

Orthogoanl states: i, j, i', j' = 1, 2 ... n, 

Physical states: cc, /?, a' , (3' — a,b... (n total) . 

The eigenstates B j can be written as linear combination of B^s: 

B i = J2 r ia B a . ( 2 . 2 ) 
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For the system composed of B° and B°, we have n — 2: 


B 1 = B\ 


Bo = B c 


and the physical states are usually written as 


{B a = pB° + qB° 

1 B b = p'B° - q'B° 


(2.3) 


or solving for B° and B°, 


\b° = c(q'B a + qB b ) 
{B° = c (p'B 0 - pB b ) ’ 


with 


def 1 

p’q + pq’ 


namely, 


r la = cq 1 , r lb = cq , 

( 2 . 4 ) 

r 2 a = cp', r 2b = -cp , 

Returning to the general case of n orthogonal states, the orthonormality of B^s can be 
expressed in terms of the physical eigenstates B a as 

Sij = {Bi\Bj) = J2 r L r jp ( B a\Bq) , ( 2 . 5 ) 

a/3 

where we have used (|2.2|) . The decay amplitude of a pure Bi state at t — 0 decaying to a 
final state / at time t is, from (|2.2| ) and (O). 

-d-Bj—»/(^) ^ ) ^iaQqf6 a (t) , 

a 


where a a f is the amplitude of B a decaying to /: 


a af = Amp(B a -> /) . 


The normalization is such that |a Q /| 2 is the partial decay rate of B a to /: 
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( 2 . 6 ) 


'y |®a/|^ 'T a • 

f 

Namely, the density of the final states, more precisely the square root of it, is absorbed into 
the definition of the amplitude. 

The time dependent decay amplitudes AB^fit) satisfy the following orthonormality 
relation [lOj. where the ‘inner product’ of A Bi ^f and Ab^j is defined by integration of 
A* B .^ f A Bj ^f over time followed by summation over all possible final states: 


roo 

E / dt A *B^f{t)A B ^f(t) 

f “ ^ 

roo 

= E Ob/3 E a a/ a ftf / dte* a (t)ef}(t) 

_ a r j 0 


a/3 / 


E/ a* af agf 


(2.7) 


— 8ij . 

In deriving the above, where we have used the generalized Bcll-Steinberger relation |8j given 


by 


J2f a af a /3f 


7a ^ 7/3 - i(m a - mp) 


— (B a \Bp) , 


( 2 . 8 ) 


together with the orthonormality of B^s fl2.5|) . Note that the relation (|2.8| ) reduces to the 


amplitude normalization condition Q2.6| ) for a = (5. While the Bell-Steinberger relation 
can be derived by requiring that unitarity is satisfied 0, it can also be derived from the 
old-fashioned perturbation theory to the lowest non-trivial order ||. The probability that a 
pure Bi at t = 0 decays to a final state / at time t is simply the square of the time-dependent 
amplitude: 

= \A Bi ^f{t)\- . 


The relation (ET 71) with i — j shows that this decay distribution conserves probability: 


E / dtV Bi ^f(t) - 1 . 


/ 


Now, take a general coherent two-body state at t — 0 given by 
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'&{t = o) = J2 c ijBiBj , 
ij 

where c tJ are arbitrary complex coefficients with 


E 



i. 


For simplicity, we will hereafter label the left and right sides of the particle pair as north 
(N) and south (S). The specific names to distinguish the two sides are irrelevant; we just 
need some labels for the two orthogonal spaces. The probability that north side decays to a 
final state In at time and the south side to a final state fs at time t$ is then 


W-/ n / s (1aVs) - \ ^2cijA B ^ N (t N )A Bj ^j 3 (ts) 


(2.9) 




From the orthonormality of the decay amplitude ( |2.7| ), one sees that this double-time decay 
distribution also conserves probability; namely, when integrated over the two decay times 
and summed over all possible final states, it becomes unity: 


r oo roc 

E / dt N / dts Tq,^f N f s (t N ,t S ) = 1. (2.10) 

r r J 0 J 0 

IN JS 

We now define the inclusive decay distribution of T to a final state /, where / can come 
from either side of the decay: 


J r _ roc 

F*—/(£) = E / dt N Tq,^f N f(t N ,t) 

In jo 
ro o 

+ E / dt s r^ff S (t,t s ), 

r JO 


which, due to fl2.10|) , satisfies 


( 2 . 11 ) 


r oo 

E / dtr^fit) = 2. 

/ 0 

The number 2 comes from the fact that the final state / can come from either side of the 
decay. The question is under what condition this is equal to the naive incoherent sum 

riSjm 4'£|ci,| 2 (r Bl ~/(t) + r B ^ f (t)) 

ij 

= ( \ C ij\ T \ C ji\ ) Ffii—j 

ij 
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which is the generalization of (|1.3|) . Using the expression of the double-time distribution 
( j2.9| ) and the orthonormality of the decay amplitude ( [2.7| ), the inclusive decay distribution 
( j2.11| ) becomes 


r *->/(*) = 


n'j 


_ -pnaive { j \ 

i^i' j 


The necessary and sufficient condition for this to be equal to is then 

GW = 0 . 




with 


( 2 . 12 ) 


Gw = XX c ij c y? + c i« c i*') 


(2.13) 


A sufficient condition for (|2.12|) to be satisfied is clearly 


Gat = 0 (for all i j- i'). 


(2.14) 


The matrix Ga> is ‘hcrmitiaiT in the sense that 

S~1 _ /^i* 

' jr M / i 

which guarantees that T^_>/(i) is a real quantity. Note also that the norm of G tl i is re-phase 
invariant; namely, when the phase of B^s are re-defined, G l% > simply changes its phase: 


B; B, 


GiV Grr> . 


Thus, the condition (|2.14|) is re-phase invariant. 

In the case of n = 2, the condition ( |2.14|) becomes Eq. (|1,4|) with G = Gu. In this case, 
it is straightforward to show that the condition G = 0 is the necessary as well as sufficient 
condition as long as y a ^ 7 &, m a ^ nib, an d coefficients p,q,p',q' are all non-zero. We still 
require that T^ a lJj(f) is correct independent of decay amplitudes a a f. The proof for general 
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case is given in the appendix. The derivation is simple if CPT is conserved in the mixing; 
namely, p' = p and q' = q. Then the condition (|2.12|) becomes 


3 ^{Gpq*) (e lat \a af \ 2 - e lbt \a bf \ 2 ) 

+2 %Gpq*) A ( e -^+- iSm ^a* af a bf ) = 0 , 

with 


(2.15) 


7a + lb x 

7+ = — 2 —’ = m « - m b . 

The three terms in Q2.15| ) have different time dependences, and thus each term should be 
separately zero. When y a ^ 7 &, this condition is equivalent to 3 %(Gpq*) = 5s(Gpq*) = 0, or 
simply 


Gpq* = 0 (CPT) . 

Thus, if both p and q are non-zero (i.e. there is a mixing), G must be zero in order for the 
naive incoherent sum to be correct. 

Let’s briefly appreciate the meaning of the condition G — 0. This is satisfied by any 
coherent state of B°B° and B°B° since then we have a = d = 0. It includes the T(4S') case 
given by (|1 . 1| ) , or a B°B° state with an even orbital angular momentum. Also, the condition 
is satisfied if b + c = 0 regardless of the values of a and d. However, the condition is not 
satisfied, for example, by the symmetric state 

B°B° + B°B° + B°B° + B°B° . 


Such a state cannot be readily produced in practice, but in principle it is possible if there 
exists an interaction with A B = 2, such as the hypothetical superweak interaction. 

To summarize, we have studied inclusive decay time distributions of coherent two-body 
states. We find that there is a set of orthonormality relations among decay time distributions 
of states that are pure and orthogonal to each other at t — 0. Using this, we have shown 
that the naive incoherent sum of single particle decay time distributions does not always 
give the correct inclusive distribution, and extracted conditions for it to be the case. Such 






incoherent sum was found to be correct for any B°B° state regardless of the orbital angular 
momentum. 
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APPENDIX: 


For n — 2, the condition (|2.12|) reduces to 

K(G'g'V) e-^Kfl 2 - ^(Gq*p) e-^\a bf \ 2 
+3? \(G*p'*q - Gq’*p ) e -(7+-«n*)t a * /a6/ 1 = 0 . 


(Al) 


Since the three terms have different time dependences, each term should separately be zero. 
The decay amplitudes are in general non-zero (we are requiring that the naive inclusive 
distribution be correct independent of decay amplitudes); thus, the above condition leads to 


K(Gg'V) =0 

< M(Gq*p) =0 , (A2) 

G*p'*q - Gq'*p = 0 

where the last condition is due to the fact that the term e l5mt samples all possible phases. 
Now we define 


G = \G\e i9 , 


def in 

s = qe 9 , 


, def 
S = 


q'e i9 . 


(A3) 


Substituting this in (|A4| ), and dividing each equation by |s| 2 , one obtains 

\g\x(G) =0 

< =0 ; (A4) 
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namely, 


|G| = 0 , or 

Thus, if \G\ 7 ^ 0,this leads to 


P P . . , 

—, - : pure imaginary, and 
s s 


P 


P 

q 


(AS) 


which means that B a and Bi, are same physical states and it contradicts our assumption 
that they have different decay rates and masses. Thus, G = 0 results from (^1). On the 
other hand, the condition G = 0 trivially leads to (^1); thus, G = 0 is the necessary and 
sufficient condition for the naive inclusive distribution to be correct assuming that y a ^ 7 &, 
Tn a 7 ^ nib and p , q,p', q' are non-zero (namely, B° and B° mix). 
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